Origin of ferroelectricity in high Tc magnetic ferroelectric CuO 
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"Magnetic ferroelectric" has been found in a wide range of spiral magnets. However, these ma- 
terials all suffer from low critical temperatures, which are usually below 40 K, due to strong spin 
frustration. Recently, CuO has been found to be multiferroic at much higher ordering tempera- 
ture (~ 230K). To clarify the origin of the high ordering temperature in CuO, we investigate the 
structural, electronic and magnetic properties of CuO via first-principles methods. We find that 
CuO has very special nearly commensurate spiral magnetic structure, which is stabilized via the 
Dzyaloshinskii-Moriya interaction. The spin frustration in CuO is relatively weak, which is one of 
the main reasons that the compound have high ordering temperature. We propose that high Tc 
magnetic ferroelectric materials can be found in double sublattices of magnetic structures similar to 
that of CuO. 

PACS numbers: 75.85.-|-t, 71.20.-b, 75.25.-j 



C/2 



cd 



o 

> 

(N 
(N 

l> 

O 

o 



X: 



Magnetic ferroelectric materials in which ferroelectric- 
ity is induced by magnetic ordering, have attracted in- 
tensive interests [l|, [^ . The strong magnetoelectric (ME) 
coupling in these materials opens up a new path to the 
design of multifunctional devices that allow the control 
of charges by the application of magnetic fields or spins 
by applying voltages. So far, almost all magnetic fer- 
roelectric materials are strongly frustrated magnets [3|- 
Frustrated magnets have very low ordering temperatures 
(^ 30 - 40 K), several times smaller than the tempera- 
tures expected from their spin interaction strengths. Low 
critical temperature is one of the major factors that limit 
the applications of these important materials. Therefore 
a new mechanism that allows high temperature magnetic 
ferroelectric materials is critical. 

Recently, CuO was found to be multiferroic at Tc=230 
K, which is much higher than the critical temperatures 
of all other magnetic ferroelectric materials [3j|. However, 
the mechanism for the high ordering temperature were 
not clear. So far, CuO is the only binary compound that 
has been found to be multiferroic Q . CuO undergoes two 
successsive magnetic phase transitions upon cooling from 
room temperature to near zero temperature. Neutron 
scattering experiments ['^ show that below rjvi=213 K, 
the spin structure is collinear antiferromagnetic (AFMl) 
[see Fig. [Ha)]. Between Tjvi and T7V2=230 K, the spin 
structure becomes non-coUinear and slightly incommen- 
surate (AFM2) [see Fig. [Hb)], with a modulation vector 
of Q= (0.006, 0, 0.017). Remarkably, an electric polar- 
ization of 160 /xC-m~^, which can be reversed by apply- 
ing an electric field of about 55 kV/m, develops in the 
AFM2 phase. The electricpolarization was attributed to 
the spiral spin structure 0, Q which was assumed to re- 
sult from spin frustration and whereas the high ordering 
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FIG. 1: Schematic sketch of the magnetic structures of (a) 
the collinear AFMl phase, (b) the noncoUinear AFM2 phase. 
The black arrows, yellow crosses, and blue circles denote the 
spin directions associated with Cu ions. The black circle in 
(a) indicates an inversion center, (c) A sketch of superex- 
change interactions Ji to J7. The single lines, double lines, 
and dashed lines represent the three types of exchange inter- 
actions between Cu ions. 



temperature is believed to come from the strong exchange 
interactions Q. 

To clarify the mechanism behind its high ordering tem- 
perature and the origin of its ferroelectricity, we carry 
out first-principles studies of the multiferroism of CuO. 
We find that CuO has very special nearly commensu- 
rate spiral magnetic structure, which is stabilized via the 



TABLE I: Calculated lattice parameters of CuO compared 
with the experimental results obtained at room temperature. 
The theoretical result is calculated using collinear antifer- 
romagnetic (AFMl) spin structure. The calculated crystal 
structure is symmertized accroding to the C2/c symmetry. 



TABLE IL Comparison of the total energies (in eV) of phases 
AFMl and AFM2, with/without spin-orbit coupling (soc) and 
with/without structure relaxation (relax). 



lattice constant 


Expt. (Refs.[18]) LSDA+U 


a (A) 


4.6837 4.5914 


&(A) 


3.4226 3.3277 


c(A) 


5.1288 5.0268 


0(6) 


0.4162 0.4110 


/? 


99.54° 100.025° 



Dzyaloshinskii-Moriya interaction [y, |7[ • The spin frus- 
tration in CuO is relatively weak, which is responsible for 
its high ordering temperature. The work suggest that 
high Tc magnetic ferroelectric materials can be found 
in materials which have double sublattices of magnetic 
structures similar to that of CuO. 

The crystal structure of CuO is monoclinic containing 
four chemical formulas per unit cell. The AFMl spin 
structure is composed of two antiferromagnetic (AFM) 
spin sublattices, in which Cu ions have the same h val- 
ues in each sublattice. The spin chains along the [101] 
direction are antiferromagnetic and are labeled chain I 
and chain II for the two sublattices, whereas the chains 
along the [101] direction are ferromagnetic. In the AFMl 
phase, all spins are aligned in the h direction, whereas 
in the AFM2 phase, chain II rotates perpendicularly to 
chain I. To accommodate the spin structures, we use a 
2x1x2 CuO supercell, which contains 32 atoms. For the 
AFM2 structure, we neglect the small incommensurate 
component of the spin structure [i.e., set Q= (0, 0, 0)] 
and rotate the spin directions of chain II 90°, so that it 
lies in the ac plane. In the AFM2 phase, the spins form 
cycloidal spirals along the a and c axis. The incommen- 
surate component of the magnetic modulation vector Q 
is extremely small, and it should not affect the calculated 
electric polarization because P oc Sj x S^- [5| . 

We perform ab initio calculations on CuO, based on 
the density functional theory (DFT) within the spin- 
polarized local density approximation (LSDA) imple- 
mented in the Vienna ab initio simulations package 
(VASP) [Sl, Isf . The on-site Coulomb interactions U are 
included for Cu ions in the rotationally invariant scheme 
introduced by Liechtenstein et al. 110!]. The spin-orbit 
coupling is taken into account in the calculation unless 
otherwise noticed. We have tried several U values. We 
present here mainly the results for [7=7.5 eV and J=0.98 
eV, which are typical values for Cu llll. Projector- 
augmented-wave (PAW) pseudopotentials [1^ with a 500 
eV plane- wave cutoff are used. A 2x4x2 /c-points mesh 
converges the results very well. We relax the structure 
until the changes of total energy in the self-consistent cal- 
culations are less than 10~^ eV and the remaining forces 
are less than 1 meV/A. 



AFMl 


AFM2 


no soc -1- relax 0.58643 


0.59588 


no soc -1- no relax 0.59481 


0.59588 


soc -f relax 


0.01083 


soc -1- no relax 0.00822 


0.01087 



To determine the crystal structure under the AFMl 
spin configuration, we relax the structure starting from 
experimentally determined structures. The room tem- 
perature crystal structure is monoclinic and of space 
group C2/c, with inversion symmetry. However, the 
magnetic structure of the AFMl phase has only P2i/c 
symmetry. Therefore, after relaxation, the crystal struc- 
ture is also reduced to P2i/c symmetry, because of the 
"exchange striction" effects. To get an idea of the ampli- 
tude of the lattice distortion, we symmetrize the crystal 
structure with P2i/c symmetry to a crystal structure 
with C2/c symmetry, following Ref.|13|. We find that 
the Cu and O ions deviate from their high symmetry 
sites by about 10^"^ A. This distorted structure preserves 
the inversion symmetry; therefore, it has no net polar- 
ization. The inversion center is shown in Fig. [Ija). An 
inversion operation about the inversion center changes 
spin chain I to chain II. The obtained structural param- 
eters are shown in Table U and the calculated structural 
parameters are in good agreement with the experiments. 

We calculate the density of states and the band gap of 
CuO of the AFMl phase. If no on-site Coulomb inter- 
action is presented for Cu, the system is metallic. For 
U=7.5 eV, the calculated band gap is 2.01 eV, which is 
close to experimental value that ranges from 1.35 to 1.6 
eV lJ-[l6J. The calculated magnetic moments for the 
Cu ion is 0.70fiB, where /i^ is the Bohr magneton, also 
in g ood agreement with experimental value 0.65-0.68/is 
171 Ha ] . The oxygen ions also have significant induced 
magnetic moments for about 0.142 fiB- 

To get the crystal structure under the AFM2 spin 
consfiguration, we constrain the spin orientations to the 
AFM2 configuration and relax the structure again. The 
relaxation also starts from the experimental high symme- 
try crystal structure. We fix the lattice vectors a, b and 
c to those obtained from the AFMl structure. We com- 
pare the total energies of the AFMl and AFM2 phases, 
shown in Table |lll After relaxation, the total energy of 
the AFM2 phase is about 0.33 meV per atom higher than 
that of the AFMl phase, which is consistent with experi- 
ments that AFM2 phase appears at a higher temperature 
than the AFMl phase [J, ll9|,|20|. Before relaxation, the 
total energy of the AFM2 phase is only about 0.08 meV 
per atom higher than the AFMl phase, mostly due to 
spin anisotropy energy. This can be seen from that if 
spin-orbit coupling is turned off, the energy difference 



TABLE III: Nonequvalent atomic positions of the AFM2 
phase in the 2x1x2 CuO super cell. 5a, Sb, and Sc are the 
deviation of atomic positions from high symmetry (C2/c) po- 
sitions. 







AFM2 




distortion (10"") 


atom 


a 


b 


c 


5a 


5b 


5c 


Cui 


0.124988 


0.249983 


0.000007 


-12.0 


-17.2 


6.9 


CU2 


0.374985 


0.749978 


0.000005 


-15.0 


-21.9 


5.4 


CU3 


0.375012 


0.249983 


0.249993 


12.0 


-17.2 


-6.9 


CU4 


0.125015 


0.749978 


0.249995 


15.0 


-21.9 


-5.4 


Oi 





0.410727 


0.125 





16.5 





O2 





0.589308 


0.375 





17.9 





O3 


0.25 


0.910733 


0.125 





22.7 





O4 


0.25 


0.089311 


0.375 





21.2 






between the two phases reduces to 0.03 meV per atom. 
This indicates that the spiral configuration is not stable 
without spin-orbit interactions, i.e., the system is almost 
invariant under a general common rotation. The remain- 
ing 0.25 meV per atom is due to the fact that the ionic 
distortion in the AFMl phase is much larger than that of 
the AFM2 phase, which will be shown below. The band 
gap of the AFM2 phase is slightly smaller than that of the 
AFMl phase. For U=7.5 eV, the band gap of the AFM2 
phase is about 1.78 eV. The local magnetic moment of 
the Cu and the oxygen ions are also slightly smaller than 
those of the AFMl phase, being 0.65 fiB and 0.107 /is, 
respectively. 

A symmetry analysis shows that the AFM2 spin struc- 
ture does NOT have inversion symmetry due to the rota- 
tion of spin chain II, which can be easily seen in Fig. [ijb) . 
The crystal structure is distorted by the Dzyaloshinskii- 
Moriya (DM) interaction [1, 0, HH, which breaks the 
inversion symmetry. The structural parameters of the 
AFM2 phase are shown in Table [1111 Compared with the 
high-symmetry structure, all the oxygen ions are shifted 
in the -l-b direction for about 7xl0~^ A, whereas all 
the Cu ions are shifted in the -b direction by a simi- 
lar amount. The ionic distortion in the AFM2 phase is 
approximately 2 orders of magnitudes smaller than that 
driven by the "exchange striction" effects in the AFMl 
phase. 

Next, we calculate the electric polarization using the 
Berry-phase theory of polarization 22| . The calculated 
total polarization is about 90 /iC m~^ in the -b direction, 
which is somewhat smaller than the experimental value 
of 160 fiC- ra~^ along the b axis. The agreement be- 
tween theoretical calculations and experimental values is 
reasonable, given that the current functionals are not ad- 
equate to treat the subtle correlation effects in magnetic 
ferroelectric materials J23| . 

In magnetic ferroelectric materials, the electric polar- 
ization can be either purely electronic or ionic. To sepa- 
rate the two contributions, we calculate the electric po- 
larization in the AFM2 spin configuration, using the sym- 



metrized crystal structure with C2/c symmetry. We ob- 
tain an electric polarization about 38 fj,C ■ m~^ along 
the 6-axis, which is the pure electronic contribution to 
the polarization. The ionic contribution is then the re- 
maining -128 /j,Cm^^. Therefore, the pure electronic con- 
tribution and ionic contribution are of the same order of 
magnitude, but of opposite sign. We also calculate the 
electric polarization using different on-site Coulomb U 
values for the Cu ions. We find that the polarization 
is very sensitive to the U values, because the electronic 
contribution and ionic contribution to the polarization 
have opposite signs. For example, if U=b eV is used, 
the polarization is reduced to about 30 /xC • m~^. We 
note that the Coulomb U is empirically chosen, which is 
somewhat unsatisfactory. However, this is the best cal- 
culations one can do presently. In the future, parameter 
free, and yet very computationally demanding methods, 
such as hybrid functional DFT [2J] can be used to solve 
this problem. 

We repeat the above calculations without including 
spin-orbital coupling. We find that the lattice distortion 
is nearly absent after turning off the spin-orbit interac- 
tion for the AFM2 phase. At the same time the electric 
polarization including both the purely electric and the 
ionic parts is nearly eliminated. This confirms that the 
spin-orbit interaction is essential to the lattice distortion 
and electric polarization in this material. 

In magnetic ferroelectric materials, the transition tem- 
peratures are predominantly determined by the mag- 
netic exchange interactions [3|. We extract the superex- 
change interactions Js of CuO using a Heisenberg model 
H = — J2ij JijSt ■ Sj —J2i(K- Sj)^ from calculated total 
energies of the different spin configurations in the sym- 
metrized C2/c crystal structure with spin-orbit coupling. 
K is the anisotropic energy due to spin-orbit coupling. 
There are 7 Js in total, which are shown in Fig. [Ijc). 
Among them, Ji is the exchange interactions between 
nearest-neighbor Cu atoms along the [101] direction, J2 
is the interaction between nearest-neighbor Cu atoms 
along the [101] direction, J7 is the interaction between 
the nearest-neighbor spins of the same sublattice along 
the b direction, and J3 and J4 are the intersublattice ex- 
change interactions. The fitted Ji= -51 meV is in a good 
agreement with J=67±20 retrieved from neutron scat- 
tering experiments [J]. The fitted value of J2=8.6 meV, 
is only about 1/6 of | Ji|, consistent with quasi-lD model 
[ii], and the fitted value of J7=9.87 meV. The fitted in- 
tersublattice couphng values J3=4.9 meV and J4=7 meV 
are weak, because the Cu-O-Cu bond angles are close to 
90° for these two Js ^].lii AFMl and AFM2 phases, 
symmetry causes additional mutual cancelation of J3 and 
J4; therefore, the energy cost of chain II rotation is low, 
as was discussed in the previous paragraph. We also cal- 
culate the next-nearest-neighbor (NNN) interactions J5, 
Jg. We find that J5 and Je are very asymmetric. J5 is 
significant with a value of -12 meV, whereas Jg is only 



2.1 meV. These values agree well with those given in Ref. 

The intrasublattice interactions Ji, J2, and J7 basi- 
cally determine the ground state spin structure AFMl. 
The next-nearest-neighbor interaction J5 further favors 
the antiferromagnetic spin chain along the [101] direction, 
whereas Jg only adds a small frustration to this config- 
uration. The major competing interactions are those of 
the intersublattice interactions J3 and J4. The weak in- 
commcnsurateness of the spin spiral caused by frustrated 
exchange interactions J3, J4 is consistent with that in 
this material the spin competition is small. We calculate 
the ordering temperature by forcing all exchange inter- 
actions to be ferromagnetic, and get Tc= 311 K by a 
Monte Carlo simulation. This temperature is only about 
1.5 times greater than the Tc of the AFMl phase from 
simulations (In RMn205 27[, the ratio is approximately 
3 - 4.), which also indicate that the spin frustration is 
weak in CuO. The lack of strong competing interactions 
in this compound might explain the high spin-ordering 
temperature of CuO. 

To summarize, we have investigated the structural, 
electronic and magnetic properties of CuO. We show that 
CuO has very special nearly commensurate spiral mag- 
netic structure, which is stabilized by the Dzyaloshinskii- 
Moriya interaction. The spin frustration in CuO is rela- 
tively weak, which is one of the main reasons that it has 
much higher ordering temperature than other magnetic 
ferroelectric materials. We propose that high Tc mag- 
netic ferroelectric can be found in materials which have 
double sublattices of magnetic structures similar to that 
of CuO. 
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